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We give an explicit realization of the "String Axiverse" discussed in Arvanitaki et. al [I] by 
extending our previous results on moduli stabilization in M theory to include axions. We extend 
the analysis of [l] to allow for high scale inflation that leads to a moduli dominated pre-BBN 
Universe. We demonstrate that an axion which solves the strong-CP problem naturally arises and 
that both the axion decay constants and GUT scale can consistently be around 2 x 10^'' GeV with 
a much smaller fine tuning than is usually expected. Constraints on the Axiverse from cosmological 
observations, namely isocurvature perturbations and tensor modes are described. Extending work 
of Fox et. al |2], we note that the observation of tensor modes at Planck will falsify the Axiverse 
completely. Finally we note that Axiverse models whose lightest axion has mass of order 10" eV 
and with decay constants of order 5 x 10^* GeV require no (anthropic) fine-tuning, though standard 
unification at 10^® GeV is difficult to accommodate. 



I. INTRODUCTION 

The dimensionless QCD 6'-anglc is a source of CP viola- 
tion in the Standard Model, highly constrained by mea- 
surements of the electric dipole moment (EDM) of the 
neutron and ^^^Hg: \6qcd\ ^ 10~^° [3]. This is the strong 
CP problem, why is 0qcd so small?. 

An elegant solution to the problem might be provided 
by the Peccei-Quinn (PQ) mechanism [3], in which Oqcd is 
promoted to a dynamical field known as the QCD axion 
{uqcd = fa'^'^ dqcd)- The axlon is charged under an anoma- 
lous global U{1) symmetry, broken by non-perturbative 
QCD effects at the mass- scale fl"'^. The potential gen- 
erated by the QCD instantons is: 
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which is minimized at 6qcd — 0, solving the strong CP 
problem. Various properties of the QCD axion have been 
considered in [5], [6]. 

However, in order to really solve the problem, this PQ 
symmetry must be of incredibly high quality. In partic- 
ular, all contributions to the axion potential from other 
PQ breaking sources must be at least ten orders of magni- 
tude suppressed compared to that from QCD instantons. 



Phenomenological approaches in field theory trying to re- 
alize such a PQ symmetry have been pursued [7]; see [8] 
for a recent approach. However, since global symmetries 
are believed to be always broken by quantum gravita- 
tional effects, it is worthwhile exploring axions in string 
theory and under what circumstances they could solve 
the strong-CP problem. 

Pseudoscalar fields with axion-like properties generi- 
cally arise in string/Af theory as the zero modes of an- 
tisymmetric tensor fields along the extra dimensions [H]. 
The number of axions is determined essentially by the 
topology of the extra dimensions. This number, like the 
gauge group and number of fermion generations, can be 
viewed as a discrete, UV boundary condition and it is 
fairly common to have hundreds, if not thousands, of ax- 
ions present. There is a separate PQ symmetry for each 
axion field, inherited from the gauge symmetry associ- 
ated with the higher dimensional tensor field. Often these 
axions pair up with geometric moduli fields to form com- 
plex chiral mutliplets in A/" = 1,D = 4 supergravity, the 
effective theory describing many string compactifications. 
Although a plethora of axions occur in string/M theory, 
they do not automatically satisfy the criteria required by 
the PQ mechanism for the QCD axion. The absence of 
PQ breaking operators is usually guaranteed to all or- 
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ders in perturbation theory because the PQ symmetries 
are remnants of (higher dimensional) gauge symmetry, 
surviving as global symmetries below the Kaluza-Klein 
scale. As in ordinary QCD, this global symmetry can be 
broken by non-perturbative effects, such as those aris- 
ing from instantons of various sorts - worldsheet, brane, 
gauge and gravitational. One has to ensure that in ad- 
dition to being valid to all orders in perturbation theory, 
the contributions to the QCD axion potential from these 
non-perturbative effects are negligible compared to that 
from QCD instantons. Axions in various string theoretic 
contexts have been studied in [TO] . 

A qualitative picture of the mass spectrum of axions in 
string theory can be given. Since the potential of a given 
axion a is generated by non-perturbative effects, its mass- 
squared will be of order j2 - where A^ ~ M^e~''^^"^\ 
where Sinst is the action of the instanton, which domi- 
nates the potential of a and 6 is a number, which charac- 
terises the " charge of the instanton" that generated the 
potential [TT] . M is given by the geometric mean of the 
supersymmetry breaking scale {V~F) and the 'fundamen- 
tal' scale, eg. the Planck scale. We will derive a precise 
formula later with the above general form. Sinst is usu- 
ally the volume of a submanifold of the extra dimensions 
in fundamental units. We will later see that fa is typi- 
cally of order the GUT scale. Therefore, the spectrum of 
axion masses will typically be exponentially hierarchical. 
Roughly speaking, with a large number of axions, one ex- 
pects a spectrum, which is uniformly distributed on a log 
scale, rather like the Yukawa couplings in M theory [12]. 
The axion masses are thus expected to span many orders 
of magnitude. These sorts of observations inspired the 
authors of [1] to propose a variety of astrophysical tests, 
which could probe the large range of axion masses from 
10~^'^eV to 10^ eV. These tests could provide evidence 
for the existence of an Axiverse. These include observa- 
tions of the polarization of the CMB, suppression of the 
matter power spectrum at small scales and the spectrum 
of gravitational waves from rapidly rotating astrophysical 
black holes. The purpose of this paper is not to discuss 
these phenomena and tests in detail, but to address in 
more depth the nature of the mass spectrum of axions 
in string/M theory and the solution to the strong CP 
problem. At the end, however, we will summarize some 
of the important astrophysical observables, which could 
either completely falsify the framework or provide strong 
evidence for it. 

The requirement of stabilizing moduli with a suffi- 
ciently large mass in realistic string compactifications 
complicates the above picture for generating masses for 
the axions. For example, simple mechanisms of stabi- 
lizing moduli with fluxes in a supersymmetric minimum 
can completely break the PQ symmetry and give axions 
masses comparable to that of the moduli [13]. Stabiliz- 
ing certain moduli by non-perturbative effects in a su- 
persymmetric minimum, such as in the KKLT idea |14j . 
also gives axion masses comparable to that of moduli; 
hence these axions cannot solve the strong CP-problem. 



In fact, it can be shown that there does not exist any 
supersymmetric minimum within Af = 1,D = 4 super- 
gravity with phenomenologically allowed values of stabi- 
lized moduli but unfixed axions (down to the QCD scale) , 
which could solve the strong CP problem 15j. In partic- 
ular, for supersymmetric vacua it was shown that either 
all moduli appear in the superpotential in which case the 
axions are very massive, or some moduli do not appear in 
the superpotential in which case the potential is tachy- 
onic [15] . 

Therefore, in order to look for QCD axion candidates 
in string theory, one should consider compactifications 
in which moduli are stabilized in a non- supersymmetric 
minimum. Such vacua have been studied in detail in the 
context of M theory compactifications on G2 manifolds 
without fiuxes [l6] [17] . In these compactifications all 
the axions pair up with geometric moduli fields in such a 
way that all the moduli superfields enjoy a PQ symme- 
try, implying that the entire moduli superpotential can 
only arise from non-perturbative effects. In a series of 
papers [161 117| it has been proven that strong dynamics 
in the hidden sector can a) generate a potential, which 
stabilizes all the moduli fields and b) generates a hierar- 
chically small supersymmetry breaking scale. Moreover, 
the potential stabilizes just one out of the many axion 
fields. Happily, the QCD axion belongs to the set of 
unfixed combinations. Higher order non-perturbative ef- 
fects, which are generically present, but not considered 
in |16l 117] because they are sub-dominant, will then sta- 
bilize the remaining axions (including the would-be QCD 
axion) at an exponentially smaller scale than the moduli 
mass scale (which is set by the gravitino mass 7713/2). In 
other words, one considers the moduli to be stabilized 
at the supersymmetry breaking scale (7713/2) with sub- 
leading non-perturbative effects being responsible for sta- 
bilizing most of the axions. Being non-perturbative, the 
resulting spectrum of axions will be distributed (roughly) 
uniformly on a logarithmic scale - as we indeed verify in 
a detailed model - thereby giving a detailed realization 
of the Axiverse. It can be shown that certain classes of 
Type IIB compactifications inspired by the moduli sta- 
bilization mechanism in AI theory [16l [17] also exhibit 
the above properties [18 . Therefore, the analysis and re- 
sults of this paper can be directly applied to the those Type 
IIB compactifications as well. For concreteness, however, 
we will study the effective theory arising from M theory 
compactifications here. 

The axion decay constants fa are constrained by vari- 
ous observations. The decay constant for the QCD axion 
j-qcd jj^ygi; |-,g bounded from below by ~ 10^ GeV because 
of processes like axion emission from stars and super- 
novae. An upper bound on all axions arises from the re- 
quirement of 'not overdosing the Universe'. The thermal 
population of light axions (< 1 eV for fa ^ 10^ GeV) 
is typically quite negligible. However, being extremely 
weakly coupled to matter, axion fields begin coherent os- 
cillations, which begin when the Hubble scale H rria, 
and the constraint arises from not storing too much en- 
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ergy density in the axion fields. This energy density is 
quadratic in the initial vev (called mis-alignment angle) 
of the axion field. The value of the upper bound ob- 
tained clearly depends on the detailed cosmological evo- 
lution i.e. matter vs. radiation domination and entropy 
releases during the oscillations. Assuming a cosmology 
characterized by a radiation-dominated phase after infla- 
tion, one gets J^^'' < lO^'^ GeV (the bound for general 
axions is given in section IV.) Thus, in order to satisfy 
the upper bound on f^'^'^ without fine-tuning the intial 
misalignement angle 9'^'^ , one has to either invoke a very 
large cycle on which the QCD gauge group is localized, 
and/or have a sufficiently low string/llD scale, imply- 
ing a large overall volume of the extra dimensions V) 
[TTj . Having a very large QCD cycle to make fa re- 
spect the above upper bound makes the strong gauge 
coupling aqcD far too weakly coupled at the compacti- 
fication scale and hence does not seem very natural. So 
the only reasonable option without fine-tuning the mis- 
alignment angle, is to have a sufficiently small string/1 ID 
scale Ms /Mil <C Mqut with a very large overall volume 
V but a small (still within the supergravity regime) QCD 
cycle such that oiqcd is reasonable [TT]. However, this 
implies that the standard unification scale ~ 10^^ GeV 
is above the string scale. Also, a low fundamental scale 
will likely lead to conflicts with current bounds on the 
proton lifetime. Of course, it is possible that the ini- 
tial misalignment angle is small (<C 1) for some reason, 
which can make a large uniflcation scale (~ 10"'^^ GeV) 
compatible with the observational bound on /o, assuming 
standard thermal cosmology. In fact, it has been argued 
that a small misalignment angle consistent with observa- 
tions could be anthropically selected [12], [131 • The M 
theory models, as we will see, (and as is assumed in [1|) 
have fa ~ Mgut lO^^ GeV. 

One of the key points of [20] is that, if the Hubble 
scale after inflation Hi > TO3/2 ^ 50 TeV then the energy 
density of the Universe is dominated by moduli oscilla- 
tions until just before BBN. This is a different cosmo- 
logical history than the one that yields the upper limit 
fa < lO^^GeV. Modufi with masses of order 50 TeV 
will decay just before BBN and, as pointed out in [5], 
the entropy released by the moduli decay can allow a 
larger /a < 5 x 10^* GeV without fine tuning. Thus, the 
anthropic fine tuning of the initial axion mis-alignment 
angle is much smaller since the upper bound is much 
closer to the GUT scale. [T] do not consider such mod- 
uli dominated cosmological histories (which are presum- 
ably generic in string theory) and in fact take the Hubble 
scale to be rather small after inflation Hj < O.lGeV. It is 
somewhat satisfying that the flne-tuning of the misalign- 
ment angles is greatly reduced in the (presumably) more 
generic case of moduli dominated, non-thermal cosmol- 
ogy. We thus extend the results of [5] to the Axiverse 
and extend the Axiverse to the presumably more generic 
case of Hj > 7713/2. We find that the non-thermal mod- 
uli dominated pre-BBN Universe is much less fine tuned 
than a standard radiation dominated pre-BBN Universe 



arising in low scale inflation. Extending the results of [2], 
which only considers the QCD axion, to the Axiverse, 
we emphasise that the observations of tensor modes by 
Planck, when combined with the existing bounds on ax- 
ion induced isocurvature perturbations, would rule out 
the string Axiverse completely, thereby requiring axion 
masses to be 'large'. 

The plan of the paper is as follows. In section |IT| the 
mechanism to stabilize the axions is outlined in detail 
and the spectrum of eigenvalues and eigenstates is com- 
puted. A simple formula for the axion masses is obtained. 
In section |III[ the effects of QCD instantons are taken 
into account and the precise method of identifying the 
QCD axion candidate among all axions is explained. It 
is shown that such a candidate typically occurs within the 



M theory framework considered. Section IV is an anal 



ysis of the axion relic abundances within the 'thermal' 
and 'non-thermal' cosmological histories, followed by a 
discussion of the allowed parameter space after impos- 
ing current cosmological constraints in section [Vj This 
is followed in section VI by a discussion of cosmological 
observables, which could either falsify or provide strong 
evidence for the framework, and also distinguish among 
the two cosmological histories. Technical details follow 
in Appendices |A| [B] and [Cj 



II. STABILIZATION OF AXIONS 

In this section we will describe the details of how, 
adding higher order non-perturbative corrections to the 
model considered in [THl [IT] , one can stabilize all the ax- 
ion flelds and calculate their masses. We work within the 
framework of supergravity in four dimensions with a su- 
perpotential W , which is a holomorphic function of the 
moduli superflelds z„i = Un + iSm whose real parts are 
the axion fields t^ and imaginary parts are the geometric 
moduli Sm- There is a real Kahler potential K depending 
on all of the Sm and whose form is given below. 

The superpotential we consider takes the form: 

00 

W = Ai4>l e''" + A2 e^^^ + ^ Ak e'^" , (2) 

where all the coefficients Ak^s are order one constants 
(because of the PQ symmetries). The first two terms 
come from the model considered in |161 I17j and arise 
from strong gauge dynamics in the hidden sector. They 
depend on only one combination of axions. The remain- 
ing axions enter through the higher order corrections rep- 
resented in the sum. These higher order corrections owe 
their origin to non-perturbative effects such as membrane 
instantons or gaugino condensates and are expected to be 
generically present in these compactifications. Note that 
in principle this is an infinite sum if one takes multiply 
wrapped instantons into account. In practice, though, 
one only requires considering as many independent terms 
as there are axions present. If the number of supersym- 
metric three-cycles in the compactification is sufficiently 
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large, it is easy to have the required number of indepen- 
dent terms. Even if the number of supersymmetric three- 
cycles is not large enough one expects non-BPS instan- 
tons to contribute to the Kahler potential, which would 
give rise to the same results qualitatively; hence for sim- 
plicity and concreteness we assume henceforth that there 
are enough independent terms present in the superpoten- 
tial above. 

The i^'s are integer linear combinations of the mod- 
uli superfields {Fk = X^ili ^K^i) ^^"^ 't'l is a holomor- 
phic composite field made of hidden sector matter fields. 
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61,2 = ^p^,Pi,2 G Z+ are the dual Coxeter 



numbers of the hidden sector gauge theory, whereas 63, 
64, ... = 27r/, / G Z+ - since the higher order terms are 
assumed to be generated by membrane instantons. This 
is consistent with these terms being higher order since 
Vk = Im(FK) > 1 is the volume of a three dimensional 
submanifold of the extra dimensions in lid units. 

There are N geometric axions ti = Re(zi) and one 
matter axion 9 (the phase of in (2|). It is also possible 
to consider a more general case with matter axions in the 
subdominant terms in We do not consider such a case 
below for simplicity. With the superpotential ^ and a 
generic Kahler potential of the form: 



K = -31og(V) 



(3) 



where V is the overall volume of the internal manifold in 
IID units, the scalar potential contains N approximately 
flat directions corresponding to N PQ symmetries, which 
are preserved by the first two terms in ([2]), while fixing 
the axion combination 1161: 



cos(xi - X2) = cos((6i - 62)^^1 - t + aS 



where xi = ■ t + aO, xk — ^k^k ■ t; K — 2,3, ... 

With all moduli and one combination of axions fixed 
by the supergravity potential coming from the first two 
terms, the remaining N would-be Goldstone bosons of 
the PQ symmetries are fixed by the next N largest terms 
in the potential, which contain N linearly independent 
combinations of axions. This is because the terms in the 
potential proportional to e~''^^^~^"'^'" cos(xi —xm) and 
^-b2Vi-bKVK cos(x2 - Xk) where M, K > 2 are much 
greater than in 6"''*'^*^ "''^^^ cos(xm - Xk)- Therefore, 
these terms fix not only the independent axion combina- 
tions xi ~ Xm and X2 ^ Xk but also effectively fix the 
combination xm -Xk,'as Xm -Xk = {xm - Xi) + (xi - 
X2) + (X2 - Xk) with xi - X2 fixed by Q. 



A. Axions - Mass Eigenvalues and Eigenstates 

In order to compute the axion spectrum, it is best to 
study the effective potential for light axions generated by 
integrating out the moduli and the heavy axion combi- 



nation. Doing that gives rise to the following: 

N+2 

Veff ^ Vo- m./^mli ^ Dk e"'-^- cos(xi - Xk) 

K=3 

Vfc: bKVK <bK+iVK+i, (5) 

where Dk are 0(1) positive numbers not important for 
the order of magnitude estimates, m3/2 = e^^^/^ (W) nipi 
and Vk is the stabilized volume of the cycle that gener- 
ated the corresponding term. 

After canonically normalizing the axion kinetic terms 
by taking the non-trivial Kahler metrics K^b = 
dz dzb ' ^ {sij^i} for moduli and matter fields into 
account, one finds the following estimate for the light 
axion masses (see Appendix [A|) : 
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where Ck are 0(1) model-dependent coefficients and the 

axion decay constants are given by = \f2K~k, with 

Kk^km = Ul^KniUim the diagonal Kahler metric ob- 
tained from the original Kahler metric by a unitary trans- 
formation U. In the notation above, the heavy axion 
combination stabilized with a mass comparable to the 
moduli is denoted by rhao ' fiance the index i labeling the 
light axions goes from 1 to TV. 

The results of [Ml [17] show that all geometric moduli 



Sj are stabilized in terms of only one parameter - (Vi), 
which is the stabilized volume of the three-cycle support- 
ing the hidden sector giving rise to the first two terms. 
This is the inverse of the coupling constant in the hid- 
den sector (Vi) = 1/ahid- {Vi) is itself determined from 
1 , (4) bi and 62 in (|2| [161 Hi]. Since Vk = E^i ^i- sj and 
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ume V are determined in terms of Sj , this 



implies that all three-cycle volumes Vk and the overall 
volume V are stabilized as well. In particular. 
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(7) 



where -/V^' -^i* ^'"'^ integers and hj are real numbers |17j . 

It is convenient to express the masses of the light axions 
in terms of the gravitino mass. Using the result for 7713/2 
in terms of nipi from [161 117] 
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one obtains the more widely applicable formula: 
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M, 



GUT 



, (9) 



where we have used fa^ ~ Mqut and V ^ 5000 [T7] . 
We have also checked that the above result ^ gives a 
very good estimate of the axion masses found exactly by 
numerical methods for a toy example with four axions 
(see Appendix [b]). 
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Since the integer coefScients in ([t]) generically differ 
by 0(1) for different three-cycles Vk, the masses will be 
distributed evenly on a logarithmic scale. Hence, this 
provides a concrete realization of the "string Axiverse" 
considered in fJ^ with a multitude of extremely light ax- 
ions (ifii <^ '713/2 ^ TeV) in a controlled and reliable 
manner. The implications of this will be discussed in 
more detail later in the paper. 

What can be said about the range of the axion masses 
possible within this framework? This range is determined 
by the allowed range of {Vj}; J — 3,N. For concreteness 
we assume that all the 67+2 = 27r. To estimate the range 
of allowed values of Vj, it is useful to note the result of 
[2T] where it was argued that the quantity rjj 



is bounded from above by 0(1) for realistic compactifica- 
tions. In many cases, it can be much smaller than unity, 
i.e. rij < 1. 

The upper bound rjj — 0(1) gives a lower bound on 
the three-cycle volume Vj, {Vj} 
|17j . the overall volume — ^ 
between sa 550 and 



A)'^'- From 
^^1/3 is generically stabilized 
1150 (in IID units), leading to 



15 ^ {V>}""" < 20. Using the resuh for -^^ 
in llTj: 
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with C2 an 0(1) constant, leads to the following upper 
bound on the light axions from ([9|: 



:O(0.1)(™3/2mp,)'/'e--{^^>'""' 
0(1) (10-^..l)eV, 



(11) 



the precise value of which depends on the stabilized inter- 
nal volume V of the compactification. As a conservative 
estimate, we take {Vj}™™ — 15, which implies m™'*^ ~ 1 
eV. 

The lower bound on light axion masses, although not 
robustly fixed by theory, can be constrained by phe- 
nomenological considerations. Successful gauge coupling 
unification in the MSSM suggests a unified gauge cou- 
pling ciq^j, sa 25, which means that the visible three- 
cycle volume Vvis ~ 25. From ([9|, this implies the exis- 
tence of at least one light axion with mass rha ~ 10~^^ 
eV. In any case, axions lighter than the current Hub- 
ble parameter Hq ss 10~^°yr~^ ~ 10"^'^ eV will not 
have started to oscillate yet and hence are not observ- 
able. Thus, for practical purposes, one has: 



10"^^ eV 



(12) 



We complete this subsection with a technical remark 
about computing the axion mass eigenstates. In super- 
gravity, the basis for the moduli and axion fields is in 
general such that the Kahler metric is both non-diagonal 
and not canonically normalized (in the sense of giving 
a canonical kinetic term for scalar fields). Therefore, 
one has to diagonalize the Kahler metric by a unitary 
transformation U as mentioned above and then rescale 



by / to get a canonically normalized kinetic term. How- 
ever, this is not the mass eigenstate basis as in this basis 
the mass matrix is not diagonal. Thus, one has to per- 
form a further unitary transformation U to find the mass 
eigenstates. Formally, one can thus relate the axion mass 
eigenstates '0a/; M = 1, 2, ...N -f 1 to the original axion 
fields tL = {t», 0}; i = l,2, .., N as: 
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MK . 



fKUj^^tL; X,i,M = 1,2, ..,7V + 1(13) 



where U diagonalizes the mass-squared matrix : 
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U\,[k^klUln', 



Vk. ^ pjf .(14) 



Note that the above analysis only studied non- 
perturbative effects arising from string instantons and/or 
hidden sector strong dynamics, and did not include the 
effects of QCD instanton contributions to the scalar po- 
tential. We turn to this issue in the next subsection. 



III. THE QCD AXION AND SOLUTION TO 
THE STRONG CP PROBLEM 



The M theory models under consideration are grand 
unified theories and, hence, one would expect the QCD 
axion to be the real part of the visible sector gauge kinetic 
function: 



N 



However, 9gcd 



(15) 



at low energies also includes the 



contribution from phases in the quark mass matrix. In 
M theory these phases also depend linearly on the ti = 
Ke{zi) [12], hence taking into account this extra phase 
shift merely amounts to modifying the coefficients Np^ 
to new coefficients N^^^. The QCD axion can then be 
written as a linear combination of mass eigenstates. Here 
we have neglected the QCD instanton contributions in 
computing the mass eigenstates, but we will include them 
shortly. 



N 



9qcd = 27rRe(F™) = ^ Nf%; U 



N N+1 

or, 0QCD = E E ^^"^^^ 

i=l K,L=1 



2tt 

7^ 



2TTt, 



N+1 



(16) 



l^KLlpL = E 7^ ' 
L = l JL 



using (13). The effective decay constants /l are defined 
by: 



JL Jr 



= 1 K=l 



(17) 



The effective decay constants f^ depend on the unitary 
matrices U and U, which in turn depend on the Kahler 
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potential. In order to get an idea of the typical size of , 
we do a simple statistical analysis by utilizing features of 
generic Kahlcr potentials consistent with G2 holonomy. 
The details are given in Appendix [C] It turns out that 
for large classes of generic Kahler potentials, /i, is close 
to the standard GUT scale Mgut, 0{lQ^^) GeV. This 
is consistent with expectation, since moduli stablization 
within this framework generically yields rripi > Mn > 
Mr:rrT [roi[T7] . This also justifies the assumption made in 
[T]. For concreteness we will use /l = Mqut = 2 x 10^^ 
GeV in our analysis henecforth. We also use the same 
value for the individual decay constants Jl for simplicity, 
which can be justified by the analysis in Appendix [C| 

We are now ready to include the effects of QCD instan- 
tons. The potential from QCD instantons alone is given 
by ([1]), expressing it in terms of axion mass eigenstates: 



V'gcd- A^^^(l-cos6l,cd) 



A"* 




(18) 



If we were to disregard the effects of other contributions 
to the potential, then the mass-squared matrix for axions 
has a particularly simple form: 



( 



A" 



The matrix (|19|) has one non-zero eigenvalue: 



(19) 



(20) 



N+l ^ 
L=l JL 

However, when we do include the effects of the other non- 
perturbative contributions, the axions will be stabilized 
as explained in section II A These vevs are generally not 
close to zero; therefore there is a danger that the linear 
combination Og^d is not fixed at a value compatible with 
the experimental upper bound \0qcd\ < 10~^°. This, in 
fact, happens if the masses of all light eigenstates ipK in 
^ are heavier than 



10" 



-10 



("^w+i) 



qcd 



eV 



(21) 



since then the QCD instanton contribution is not strong 
enough to shift the vevs of V'/f's away from the values set 
by the supergravity scalar potential for light axions ([s]). 
Interestingly, rrig^xp is close to the mass of the axion for a 
volume of order l/acuT- 

This also implies that, since the QCD axion is repre- 
sented by a linear combination of the mass eigenstates, in 
order to achieve compatibility with the experimental limit 
on 9 qcd, it would be enough if just one of the mass eigen- 
states {iPk} contained inside the linear combination 9qcd 
was light compared to m^xp- In this case, although the 
QCD effects will not perturb the vevs of the heavy eigen- 
states, the vev of the lightest eigenstate contained inside 



dqcd will get readjusted to a new value such that the to- 
tal linear combination 6qcd is minimized at 6 qcd < 10"^*^. 
In addition, the mass of this lightest eigenstate will be 
almost entirely determined by QCD instantons. 

From the analysis of the axion spectra in section |II A| 
we find that the axion masses are distributed linearly on 
a logarithmic scale. Thus, by choosing natural values 
of microscopic parameters in Q and ([S]), we expect 
to find among the mass eigenstates {ipx} a number of 
very light modes with masses smaller than m^xp, imply- 
ing that the experimental upper bound on 9qcd can be 
easily satisfied. In particular, taking into account the ef- 
fects of the QCD instantons in the full mass-squared ma- 
trix will drastically modify the lightest eigenvalue, which 



now gets a mass m^' 



qcd 



r, and will modify all eigen- 



states with mass less than rricxp- The eigenstate with the 
mass m'^'^ will then be uniquely identified as the QCD 
axion for all practical purposes. 

In Appendix [B] we consider a toy example in detail in 
which we compute the eigenvalues and eigenstates of the 
axions numerically. The numerical results are completely 
consistent with the general results described above. 



IV. COSMOLOGICAL EVOLUTION AND 
RELIC ABUNDANCE 



We now study the cosmological evolution of the N ax- 
ions and compute their relic abundance. However, doing 
so requires a knowledge of the cosmological history as well 
as a knowledge of the mass spectrum of moduli vis-a-vis 
the Hubble parameter during inflation Hi . Although Hi 
is not determined from observations at present, it is pos- 
sible to get an idea about its magnitude. It is known 
that the slow-roll parameter for simple models of infla- 
tion - e = TOpj I where V is the slow-roll inflaton 

potential and V is the derivative of V with respect to 
the inflaton field, can be written in terms of Hi as: 



10^ 



mpij 



(22) 



using the experimental value of primordial density per- 
turbations 5p/p ~ 10"^. This implies that Hi < 
10"^ nipi using the fact that e < 10"^ is required for 
~^60 e-foldings of inflation to solve the flatness and hori- 
zon problems. This is the standard flne-tuning required 
in slow-roll inflation models. A smaller value of Hi 
than above will make e even smaller, implying a larger 
flne-tuning for e than is required for inflation. In this 
paper, we will take an unbiased approach and discuss 
both cases with a large Hi{> Mmoduii) and with a small 
Hi{< Mmoduii), since the cosmological consequences are 
qualitatively different. We discuss cosmological observ- 
ables, which could distinguish between the two situa- 
tions. 
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A. Relic Abundance in Non-Thermal Cosmology 

We first assume that the Hubble parameter during in- 
flation is large, i.e. Hi > MmoduU = O {1713/2), with 
'^3/2 ~ 10 TeV required for low-scale supersymmetry 
with gravity mediation, which is natural within this 
framework |161ll7j. For notational convenience, this cos- 
mological scenario will be dubbed "non-thermal" cosmol- 
ogy. In this case, there generically exists at least one 
set of moduli Xught that are lighter than Hi, which in 
turn will generically be displaced from their latc-time 
minima during inflation. After the end of inflation, the 
Hubble parameter keeps decreasing and will eventually 
become comparable to the masses of these moduli, at 
which time they will start coherent oscillations. Since 
these oscillations scale like matter, they will quickly dom- 
inate the energy density of the Universe. The lightest of 
these moduli will start oscillating last. Since the mod- 
uli are very weakly coupled to the visible sector, they 
will decay long after all the moduli have started oscil- 
lating. The requirement that the decay of the lightest 
modulus Xq gives rise to a reheat temperature greater 
than a few MeV for successful BBN puts a lower bound 
on the mass of Xg. Within gravity mediation, one finds 
mxo ^ 0(10) TeV, the precise value depending 

on model-dependent details. It was shown in [20] that 
both the moduli and gravitino problems can be naturally 
solved within this framework. 

In order to compute the relic abundance of axions in 
such a framework, it is important to consider two era's: 
before and after the moduli have decayed. During the 
first period the Universe is moduli dominated. This pe- 
riod then gives way to a radiation-dominated era after the 
decay of the lightest modulus. It is this latter period in 
which BBN and later cosmological events such as recom- 
bination, matter-radiation equality, and growth of struc- 
ture take place. The computation of the relic-abundance 
of axions in these two regimes is very different, which we 
now study. 



1. Radiation- dominated Era 

If the mass of a given axion is such that it starts 
oscillating after the lightest modulus decays, i.e when 

rha < = 0{1)^^^^, the standard computation of the 

relic abundance is applicable. This gives the following 
boundary value: 

m = 0(1 - 10) X 10"^^ eV (23) 

For axions with masses below m, the axion relic abun- 
dance can be computed as: 

/ - \ 2 

/ f \ / m \ 1/2 

^,2 _ Q Qg ( Jak \ I '"'ak \ ' lal, 



2 X 1016 GcV 



10-2" eV 



mm 



where x is an 0(1) fudge factor to take into account pos- 
sible effects from anharmonicity, etc. The masses of these 



axions are due to string/membrane instantons and are 
much smaller than that of the QCD axion (« 3 x lO"!" 
eV), hence they are not expected to receive corrections 
from finite temperature effects. It is important to note 
the dependence of the relic abundance on fa^ and rha^. ■ 
The relic abundance increases by increasing the mass and 
the decay constant. The upper bound on the relic abun- 
dance therefore implies that for /^^ « 2 x lO^^ GeV and 
(Of^) = C(l), there is an upper bound on the mass of the 
axion: 



m, 



(std) 



0(1) X 10"^° eV 



(25) 



Thus, from (23) and (25), the misalignment angle {Oj ) 



has to be fine-tuned for axion masses between lO^^o < 
ma. < 10-14 eV. 



2. Moduli-dominated Era 

For axions with masses rhaf. > 10~^'^ eV, the relic 
abundance is determined by a different computation since 
the Universe is moduli-dominated. The abundance can 
be readily computed 2 : 



^ak h' 



0{l) 



0(10) 



rpXo f2 ^ 

^RHJak 

mil (3.6 eV)^ 

I fok 

\ 2 X 10i6GeV 



{ol)x 



(26) 



rjiXo 

IMeV 



^ik)x 



Note that the relic abundance is completely independent 
of the mass of the axion; apart from (Of^) it only depends 
on the (effective) decay constant, which is approximately 
Mgut = 2 X 10^6 GeV within our framework, and the fi- 
nal reheat temperature T^fj , which is more or less around 
1 MeV. Thus, the relic abundance for all axions with 
masses rhak ^ 10- eV is independent of their masses 
and only depends on the misalignment angle (O'j^)- Note- 
ice that, within our framework, the mass of the QCD 
axion (the mass eigenstate that dominantly gets its mass 
from QCD instantons) automatically lies in this region, 
mf^ « 3 X 10-1° eV > lO-i'^eV. Hence its relic abun- 
dance is just determined by {9j)aqod. In order to be con- 
sistent with the WMAP upper bound {nDMh"^ < 0.11), 
this requires: 



< 10' 



(27) 



which implies a modest fine-tuning in the intial condi- 
tions. We will discuss this in detail in section [V] 



B. Relic Abundance in 'Thermal' Cosmology 

The cosmological framework studied above can be con- 
trasted with one in which Hj < MmoduU — 0(10) 
TeV, which will be termed as "thermal" cosmology for 
notational purposes,. This corresponds to a situation 
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with a much smaller Hi and hence a much smaller re- 
heat temperature after inflation T^jj. Now, a radiation- 
dominated phase follows after reheating (from inflation) 
at T^jj and continues until matter-radiation equality at 
Teq- In this case the relic abundance of all axions is 
computed as it was done in section [IV A 1[ giving rise to 
the same equation: 



O ^ /i-^ = 0.06 



2 X 1016 GeV 



10- 



'eV 



1/2 



{elm 



For the QCD axion, finite-temperature mass effects for 
the QCD axion with. T^aoco ~ 3 x 10-1° eV will modify 



the dependence on {/afc,Was,} compared to that in (24 1 
Thus, for the QCD axion one finds [2]: 



^''o.Qcn"' 



4.5 X 10* 



far 



7/6 



2 X lOieCeV 



Note that (28) means that the relic abundance for each 



axion increases with the mass of the axion in contrast 
to that within non-thermal cosmology, giving rise to an 
upper bound for the axion mass rn'reUc 10"^'^ eV, above 
which (6*/^) has to be tuned. 



V. OTHER COSMOLOGICAL CONSTRAINTS 

In addition to the relic abundance constraint, there 
are other constraints arising from the presence of light 
axions. The first is the presence of an isocurvature com- 
ponent of temperature fluctuations (aa) and the second 
is the presence of a non-Gaussian component of tem- 
perature fluctuations, both arising from fluctuations of 
the axions during inflation. These constraints have been 
studied earlier in various contexts. For example, [2] has 
studied the constraints for the QCD axion with large /^'^'' 
for both thermal and non-thermal cosmological histories. 
A recent paper |22j studies the consequences for iscourva- 
ture fluctuations within the above assumptions for stan- 
dard cosmology. [23] has studied the constraints for the 
QCD axion with thermal cosmology as a function of the 
decay constant /^'^'^. [24] has studied constraints for the 
QCD axion for both standard and non-standard cosmo- 
logical histories assuming that the QCD axion comprises 
all of DM. [25J has studied constraints arising from mul- 
tiple axions, but with different assumptions than the one 
considered in this paper. In our analysis, we have fixed 
fa^ to be the GUT scale as it is the natural scale in the 
framework considered. We then generalize the analysis 
of [2| for the QCD axion to a situation with many axions 
roughly distributed evenly on a logarithmic scale. 

The observables above depend on the axion relic abun- 
dances {Q.ah?) and the Hubble parameter during infla- 
tion {Hi). Furthermore, the gravity wave contributions 
to temperature fluctuations also depend on Hi. It turns 
out that the bound on a non-Gaussian component in the 



CMB does not give rise to any additional constraints on 
the parameter space, so we consider the following: 



nah" = > na^e < 0.11 



N 

E 

k=l 
N 



^ ( (^Qt/^m)^ ^ 2 Inn 



(29) 



cri ) < 0.072 



Qt 



Hr 



StT TJlpl 



< 9.3 X 10" 



where we have used the latest bounds from WMAP5 [3B] . 
Note that the quantity {Oj^) appearing in the expression 
for the relic abundance of axions is given by: 

/n2 \ _ /)2 , _2 



where ag. 



Hi 

2^ fa. 



(30) 



The observed upper bounds on the relic abundance 
(fia/i^), tensor modes {Qt) and isocurvature fluctuations 
{aa) provide a constraint on the 2N + 1 microscopic 
parameter space - {rha^^Oi^., Hi; k = l,iV} in general, 
where 0/^ is the initial mean misalignment angle of the 
axion ak. More precisely, a given spectrum of axions 
{rria^} imposes constraints on the parameters 

From our understanding of the mass spectra of axions 
in section |II A[ it is possible to dramatically reduce the 
number of parameters, as follows. Using (11) and (12), 



lO^^^eV < m^, < leV; 
39 > Vj > 15; 



k 

J = 3, TV 



1, 



(31) 



Note that for non-thermal cosmology, the boundary be- 
tween the moduli-dominated and radiation-dominated 
regimes given by m = 0(1 — 10) x 10"!^ eV in (23l, 
corresponds to Vj w 25. Now, since rhak e~'^^'^+^ 
and Vk+2 varies by 0(1) for different three-cycles in the 
compactification, we assume that there are 0(1) axions 
in each e-folding between ^ 1 eV and ^ lO^'^'^ eV, cor- 
responding to 15 < Vk+2 ^ 39. This is expected to 
be true provided the total number of axions N is suffi- 
ciently large. For simplicity we also assume that the ini- 
tial mean misalignment angles of all axions are roughly 
equal: 9i^ ~ Oi^, k = 1,N. 

With these assumptions, the number of microscopic 
parameters is reduced to four for non-thermal cosmo- 
logical evolution - {9iq, Hi, Ni, N2}, where A^i and N2 
are the number of axions with masses corresponding to 
Vk+2 hi the ranges 15 < Vk+2 ^ Vk+2{~ 25) (moduli- 
dominated regime) and 25 < Vk+2 ^ 39 (radiation- 
dominated regime) respectively. For thermal cosmology, 
the number of parameters is just three - {Oi^, Hi, Nstd}, 
where N^td is the number axions with masses correspond- 
ing to Vk+2 in the range 15 < Vk+2 ^ 39. In particular, 
with one axion in each e-folding between ~ 1 eV and 
- 10-33 eV, A^i « 10, N2 « 14, and Nstd = 24. 

Figure [T] shows the effect of the cosmological con- 
straints on the two dimensional parameter space 
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Then, the number of axions Nstd in the entire mass range 
correspond to 15 < Vfc < 39, implying Nstd = 24. The 



-4 -3 -2 
LogioLfio] 



FIG. 1. Allowed microscopic parameter space (unshaded 
region) in the {dig,Hi} plane for A^i = 10, N2 = 14 
with a "non-thermal", moduli dominated cosmological his- 
tory {Hi > Mmoduii) after imposing the current bounds on 
tensor modes, isocurvature fluctuations and the overall relic 
abundance. Contours for three allowed values of the isocur- 
vature fluctuations aa are also plotted. 



{digjH]} for non-thermal cosmology with A^i — 10 and 
N2 = 14. The range of Hj is chosen to be between 10^ 
and 10^^ GeV. The lower bound is so chosen because the 
non-thermal cosmology occurs when Hj > MmoduU- It 
can be seen from the figure that the bound on tensor 
modes excludes Hj > 3 x 10"'^* GeV, while the combina- 
tion of the isocurvature and relic density bounds imply 
Hi < 10^° GeV and 9io < IQ-^. Thus, a modest fine- 
tuning (percent level) in Oi^ is required. 

It is instructive to compare the results obtained in 
Figure [T] with those in the thermal cosmology, i.e. in 
which Hi < MmoduH- As explained earlier, in this case 
the moduli are not displaced from their late-time min- 
ima during inflation, and a standard thermal history 
with a radiation-dominated phase follows after the end 
of inflation. With Mmoduii ™3/2 }^ 10 TeV as in 
our framework. Hi is bounded from above by approx- 
imately 10^ GeV. Assuming that the reheating process 
after the end of inflation is efhcient, this implies that 
the reheat temperature after inflation T^^^ ^ 2 x 10^^ 
GeV. Then, the thermal abundances of the gravitino and 
axinos/modulinos (also with mass sa 771^/2 in our frame- 
work [2^) are such that they do not overproduce LSPs 
by their decays [371 155]. Moreover, since 1713/2 > 10 TeV, 
they decay before the onset of BBN. 

Figure |2] shows the constraints on the microscopic pa- 
rameters {9ig,Hi} taking into account the above effects 
with the same mass distribution of axions, i.e. with one 
axion in each e-folding between ~ 1 eV and ~ 10~^^ eV. 
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FIG. 2. Allowed microscopic parameter space (unshaded re- 
gion) in the {^/p,///} plane for Natd ~ 24 with a "thermal" 
cosmological history (Hi < Mmoduu) after imposing the cur- 
rent bounds on tensor modes, isocurvature fluctuations and 
the overall relic abundance. Contours for three allowed values 
of the isocurvature fluctuations are also plotted. 

relic abundance bound requires 9ig < 10^^, while the 
isocurvature bounds are automatically satisfied by im- 
posing the relic abundance bound and the requirement 
Hi < MmoduH- This further implies that isocurvature 
fluctuations can only be observed within thermal cos- 
mology with aa ^ 10~^. 

The differences between Figures [T] and [2] are clear. In 
Figure [2] the relic density bound implies a fine-tuning in 
9ig , which is a little less than three orders of magnitude 
worse compared to that for Figure [l] This is because the 
relic abundance is a monotonically increasing function of 
the axion mass in Figure [2] for all axions, while it is inde- 
pendent of the axion mass (for rhaf. ^ 10"^^ eV) in Figure 
[1] This is crucial because the largest contribution to the 
relic abundance in Figure [2] comes from these heavier ax- 
ions. From (22), the requirement Hi < Mmoduii suggests 



that the fine-tuning in e (at least for simple inflation- 
ary models) is much worse than that in Figure [l] Our 
results are consistent with the earlier results of reduced 
fine-tuning for a single (QCD) axion coming from en- 
tropy production due to late decay of scalar condensates 
[2], [29], and generalizes those results to the case with 
a plethora of axions. Finally, the magnitude of allowed 
isocurvature fluctuations is about five orders of magni- 
tude smaller for thermal cosmology compared to that for 
non-thermal cosmology. 

In the above, we have assumed the existence of one 
axion in each e-folding between ~ 1 eV and ^ 10"^'^ eV. 
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However, it could happen that Vk for different three- 
cycles in the internal manifold scans less finely (but still 
varies by C(l)), giving rise to say, one axion in every 
ten e- foldings. In addition, as mentioned below (11), de- 



pending upon the details of the compactification V could 
be stabilized at values close to its upper bound, lead- 
ing to a larger lower bound for Vk and hence a smaller 
TO™°^. These effects will make {A^i,iV2} (Nstd) smaller 
than that assumed in Figure 1 (Figure 2), and could help 
in relaxing the constraints on {9ig,Hi}. 




FIG. 3. EfTect on allowed values of microscopic parameters 
{dig, Hi} by decreasing {A''i, for non-thermal cosmologi- 
cal history {Hi > M^oduu)- Left: iVi = 10, iV2 = 14; Right: 
iVi = l,iV2 =2. 



In Figures [3] and |4j we show the effects of decreasing 
{7Vi,iV2} from {10,14} to {1,2}, and N^td from 24 to 
3, respectively. Within non-thermal cosmology this has 
the effect of relaxing the constraints on {fl/p, i?/}, as seen 
from Figure [3] However, the constraints for thermal cos- 
mology shown in Figure |4] are essentially unchanged. 
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FIG. 4. EfTect on allowed values of microscopic parameters 
(S/p , Hi} by decreasing Nstd, for thermal cosmological history 
{Hi < M^oduu). Left: Nstd = 24; Right: Nstd = 3. 



The above analysis shows that a String Axiverse with 
a large Hubble parameter during inflation Hj > M^oduH 
is much less constrained than in the alternative case. 



A. Consequences 

It is important to understand and appreciate observ- 
ables which could falsify the entire approach as well as 
distinguish among the two cosmological histories. From 
the third equation in ( 29 1 and from Figures [l]and[2j it can 
be seen that the observation of tensor modes in the future 
by PLANCK requires a reasonably large Hj > 3 x 10^^ 
GeV, which is conclusively ruled out within our approach 
for both cosmological histories. Hence, an observation 
of tensor modes in the near future in general, and by 
PLANCK in particular, will rule out the String Axiverse. 
This is a similar conclusion to but now it applies to 
the entire Axiverse. Note also that isocurvature fluctua- 
tions can effectively discriminate between the "thermal" 
and "non-thermal" cosmological histories considered in 
this paper. As can be seen from Figures [T] and [2] non- 
thermal cosmology can give rise to isocurvature fluctua- 
tions just below the observed bound q-jjo"-'**'' < 0.072, 
while thermal cosmology predicts isocurvature fluctua- 
tions which are vastly suppressed, a^*'' ^ 3 x 10~^. 
Therefore, an observation of isocurvature fluctuations in 
the near future will rule out a String Axiverse with ther- 
mal cosmology . On the other hand, although a lack of 
observation of isocurvature fluctuations in the near fu- 
ture will not rule out non-thermal cosmology within the 
approach, it will disfavor it. 



Constraints from Production of Light Axions 
from Other Sources 



The relic abundance of light axions in the previous sec- 
tions has been computed assuming that the axions act as 
coherent classical fields with zero momentum. However, 
in general there are two other contributions to the axion 
relic abundance: 

• Light axions could be produced during interactions 
among particles in the thermal plasma created dur- 
ing reheating after inflation. 

• Light axions could be produced directly from the 
decay of moduli (scalar fields) with masses ~ "1.3/2- 

The production of light axions in the thermal plasma 
has been studied in 30J for "thermal" cosmology. It turns 
out that axions with fa^ ~ 2 x 10^^ GeV interact so 
weakly with the thermal plasma that a thermal popula- 
tion of axions never results. 

Within "non-thermal" cosmology, these thermal ax- 
ions, even if present, will be vastly diluted by the decay 
of moduli. So they are completely irrelevant. However, 
in this case light axions can be produced from the decay 
of coherently oscillating scalar fields displaced during in- 
flation. The axions thus produced contribute to the total 
number of effective neutrino species N^^^ — N^, + AN^, 
for which there is an upper bound from BBN due to 
"^He overproduction as well as from CMB measurements 
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(AA^^ < 1) [mUH]. Hence, it is important to check if 
these bounds can be satisfied within our framework. 

The relevant scalar field to consider is the one that 
is the lightest and decays last, since the axions possibly 
produced from the decay of heavier moduli will be vastly 
diluted by entropy production of the lighter ones. As 
noted earlier, within the scheme of moduli stabilization 
considered here, A'^ — 1 moduli are stabilized with masses 
w (1 — 2) 7713/2 EH]- In cases with matter axions present 
in the subdominant condensates, it turns out that the 
matter axions are also stabilized with masses « TO3/2- If 
it happens that the lightest scalar field Xq is an axion 
rather than a geometric modulus, then it cannot decay 
to two light axions because of the derivative nature of the 
axion coupling. In this case, there is negligible produc- 
tion of ultra-relativistic light axions. On the other hand, 
if Xq is a geometric modulus, then a tree-level decay of 
this field to two light axions is possible via the operator 
: Ck Xq d^dkd^^dk for some model-dependent coefficient 
Cfc. Since the decaying scalar field dominates the en- 
ergy density of the Universe, the yield = AN^ of light 
axions is given by: 

(32) 

where Ba is the branching ratio of Xq to all axions - 
^0 ~^ o,k cik', k — 1, N . Here Tjih is the reheat tempera- 
ture from the decay of the lightest modulus Xq, which is 
0(1) MeV for Mxg > 10 TeV. The bound is thus easily 
satisfied. 



VI. OBSERVABLES 

We finish the paper by returning to the observable con- 
sequences of the String Axiverse that were discussed in 
[1]. The basic prediction for the Axiverse spectrum we 
find is that the masses are distributed roughly evenly on 
a logarithmic scale between ~ 1 eV and ^ 10"^'^ eV. 

As explained in [T], the mass spectrum of axions can 
be divided into four windows as far as observable ef- 
fects are concerned. If there are axions in the window 
gY < < 4 X 10~^® eV, which couple appre- 

ciably to E ■ B, they could cause a rotation of the po- 
larization of the CMB. Could such axions arise within 
the framework considered here? Unfortunately not. The 
M theory models have Standard Grand Unification and 
hence the axion that predominantly couples to i? • _B is 
the QCD axion, whose mass is too large (w 3 x 10^^° 
eV). Since the QCD axion is a linear combination of 
all axion mass eigenstates one might obtain couplings of 
lighter axions to E ■ B, however these are suppressed by 
0{ "^heavy )^; coufirmcd in the toy example considered 
in Appendix [Bj Hence, one does not expect a rotation of 
the polarization of the CMB. 

Axions in the mass window 10~^* eV < < 3 x 

10~^® eV can give rise to step-like features in the matter 



power spectrum at small scales. The physics underlying 
this effect is the following: for very light scalar fields, 
there is a mass-dependent pressure term in the Euler 
equation governing the density fluctuations, which gives 
rise to a modified Jeans length, below which density fluc- 
tuations do not grow. It is very similar to the suppression 
from free-streaming due to light neutrinos. Such axions 
can naturally arise within our framework. 

Axions in the mass window 10^^* < iha^ < 10^^° cV 
can form exponentially growing bound states with rotat- 
ing black holes and hence significantly affect their dy- 
namics by graviton emission that carries aways the black 
hole's angular momentum [T. This causes a spin-down 
of the black hole, and occurs for black holes lighter than 
^ lO''^ Msoiar, (corresponding to axion masses heavier 
than ~ 10~^® eV), resulting in gaps in the mass spectrum 
of rapidly rotating black holes. A recent paper [32] stud- 
ies other interesting signals such as gravitational waves 
and gamma rays emanating from this axion-black hole 
bound state under particular conditions. These can sup- 
posedly probe axions with masses upto ~ 10~^° eV, the 
mass of the QCD axion. Since the effect only depends 
on the mass of the axion, the signal is the same for both 
cosmological histories considered. 

Finally, if axions with masses in the range 10^^*^ < 
^ofc ^ 1 eV have an appreciable coupling to E ■ B, 
they may be detected by their decays to photons in our 
galaxy or by spectroscopic effects in compact astrophysi- 
cal environments of magnetars, pulsars and quasars aris- 
ing from photon-axion conversion in the strong magnetic 
field present inside these objects |6j. The latter may be 
crucial in probing large axion decay constants fak ^ 10^ 
GeV |33j . Do we expect such axions within our frame- 
work? As explained earlier, before turning on QCD in- 
stanton effects the axionic partner of the modulus mea- 
suring the QCD gauge coupling is an 0(1) linear com- 
bination of all axion mass eigenstates. From the argu- 
ment in the paragraphs above, we saw that after turning 
on QCD effects eigenstates, which are much lighter than 
ml'^'^ ~ 10~^° eV do not couple appreciably to E ■ B. 
However, eigenstates that are much heavier than 10^^" 
eV can still couple appreciably to E ■ B since they are 
hardly affected by QCD instantons. This can also be 
checked from the explicit example studied in Appendix 
[b] Thus, we expect a few axions in the range 10~^" eV 
< TTia^ ^ 1 eV coupling appreicably to E ■ B. The feasi- 
bility of such signals, however, depends on the strength 
of axion-photon coupling (~ fa^) and the reduction of 
backgrounds. 

It is interesting to note that if the lightest axion in the 
Axiverse has a mass corresponding to a GUT instanton, 
i.e. m^j. > 10~^^ eV for all axions and if additionally the 
axion decay constants fa are all of order 10^'' GeV, then 
there is no flne-tuning required at all to satisfy the cosmo- 
logical constraints in the non-thermal, moduli dominated 
cosmology. In this case, the fundamental scale is also of 
order 10^^ GeV. Although not fine tuned, this case may 
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be difficult to reconcile with both unification and bounds 
on the proton lifetime within standard GUTs. It has 
been argued however that within orbifold GUTs it may 
be possible to satisfy both unification and proton decay 
bounds with a lower scale [34) . 

This framework generically leads to a significant com- 
ponent of dark matter in the form of axions. This is 
expected to be true for both M theory compactifications 
considered here and the classes of Type IIB compactifi- 
cations considered in [18]. For the M theory case with 
Hj > M^o^uU in particular, when combined with the re- 
sults of [IF, which show that, if stable, there is also a 
significant component of neutral Wino dark matter, we 
are led to predict two significant sources of dark matter. 
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Appendix A: Mass Eigenvalues 

In this section, we estimate the mass eigenvalues of the 
light axions. We will consider the superpotential in ^ 
with TV geometric axions and one matter axion (0i in 
the first term in ([2])). To that effect, let us define the 
following A'^ -|- 1 linear-independent axionic variables 



K 



Xi-XK+i, K = l,N+l. 



(Al) 



Recall that Ti has already been fixed by Q at = 
TT -|- 27rni, where rii e Z. By examining ([S]) it is easy to 
see that in order to minimize the potential, the remaining 
N axions Ti , i > \ will be fixed at the values T° = 27rni, 
where e Z. Then, the corresponding N x N Hessian 
matrix for the light axions evaluated at the minimum is 
approximately diagonal and positive-definite: 



Ti=T\ 
K/2 



Possible off-diagonal entries could have come from the 
terms in ([s]) proportional to 6"^'=^'="''™^'", where k,ra> 
2. However, as argued in section [Hj such terms are safe 
to neglect in the "double condensate" regime. 

Note that the "heavy axion" Ti is an eigenvector of 
the A^ + 1 X A^ + 1 Hessian Vab , which should in princi- 
ple be included along with the light axions. Clearly, the 
linear-independent combinations {Tk} are eigenvectors 



of the Hessian. However, they do not constitute an or- 
thonormal basis. Indeed, the "vectors" {Tk} are related 
to the "basis vectors" {ti, 6} by a general GL{N + 1, i?) 
linear transformation, which is generally not an orthog- 
onal transformation. To construct an orthonormal basis 
we can use the Gram-Schmidt process as follows: 

C/i = Ti , (A3) 



C/2 -Ta 



{Ui ■ T2) , 



Un+1 — Tj 



U, 



Ul 



\Ui\ 



N+1 



,U2 



N 

J = l ' 
U2 

\U2\ ' ■ 



u. 



■{Uj ■ Tn+i) , 



Un+1 = 



Un+1 



\Un+i\ ■ 

Using the orthonormal basis {uk} we can construct an 
orthogonal matrix R G 0{N + 1) whose j-th column cor- 
responds to the components of Uj in the basis {ti, 6}. 
Hence, the eigenvectors {u^} are related to {ti, 6} by the 
O(A^-l-l) rotation. Let us now compare the eigenvalues of 
the Hessian corresponding to the Tks with the eigenval- 
ues obtained in the orthonormal basis {uk}- Note that 
since the transformation between the Tr-s and the ukS is 
generally not orthogonal, the corresponding eigenvalues 
are not necessarily the same. Using ( A3 1 one can perform 



an inverse transformation to express the T^s in terms of 
the orthonormal eigenvectors {uk} 

Ti =cim, (A4) 
T2 = C2U2 + C21U1 , 

= C3U3 + C31U1 -I- C32M2 , 

n 

Tn+i = cn+iUn+1 + ''^C(^N+i)jUj ■ 
From the above, we see an important property 



dTi 



K 



du 



= 0,VAf > K 



dV dV 

= ck- 



M 



duK 



dT, 



K 



N+1 

j=k+i 



dV 
dfi 



(A5) 

Furthermore, due to the exponential hierarchy of the 
terms retained in the sum in ([S]) we generically expect 



dV dV 



dV 
duk 



Ck 



dV 
dfk 



(A6) 



(A2) Thus, we conclude that up to some multiplicative factors 
of ~ the eigenvalues of the Hessian in the or- 

thonormal basis {ui} are essentially the same as the ones 
in the {T^} basis! Since in this basis the Hessian ma- 
trix ( |A2| is approximately diagonal, all we need to do is 
determine the axion decay constants fi by finding a uni- 
tary transformation Uki, which diagonalizes the Kahler 
metric for the axions and then find the eigenvalues of the 
matrix: 

1 



fi+2f- 



'UijrS^mnUnj 



(A7) 



i+2 
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In the Kahler metric, the ofF-diagonal entries are sup- 
pressed by a factor hi ~ 0{\/N) or hi/ Si relative to the 
diagonal ones [T7] . Therefore, treating the Kahler metric 
as diagonal is justified, which will be confirmed by an ex- 
plicit numerical example in appendix [Bj The main effect 
on the eigenvalues of the Hessian then comes from the 
rescaling by the inverse decay constants. Thus, a reliable 
order of magnitude estimate for the masses of the N light 
axions is given by: 



2 _ "W2^^(X)/2, 
fi+2 



(A8) 



since fk ~ Mqut- Here Ci is a coefficient of 0(1). This 
is the expression used in (|6|. 



Appendix B: An Explicit (toy) Example 

We now study an explicit (toy) example with N = 2 
geometric moduli and a single complex matter field to 
illustrate the main features of the axion mass spectrum 
and its decay constants. The Kahler potential, superpo- 
tential and gauge kinetic functions are taken as: 

if = -31nV-H47ri/3^i^: V = 4t:^^^ s? si , 

9 / p ^- 27r f 1 ^2tv f2 ■ 2-k f 3 

+ A^e Pi ^ , 

For the following choice of parameters: 

Ai^27, ^2 = 2.27665, A3 ^ 3 , A4 ^ 5 , 
Pi =27, P2 = 30, F3=4, P4-3, 

we obtain 

si w 48.82 , S2 w 24.41 , w 53.81 , 

ti w 5, ^2 ~ -10, 6I1 w -157r. (Bl) 

The geometric moduli si, S2 and the meson 05 form three 
mass eigenstates with masses 



in section [Tl A| and in appendix [X] that the Kahler metric 
in the {t, 9) basis is essentially diagonal. From the cor- 
responding eigenvalues we now compute the axion decay 
constants Jl = = \/ 2Kl, labeled in the same order 
as the eigenvectors corresponding to the columns of U 
above: 



rripi 



(3.03xl0-^ 6.05xl0-^ 1.22). (B4) 



These decay constants are then used to rescale the 
above eigenvectors 



Ukl ^ UklIl, no sum overL , 



(B5) 



to obtain canonical kinetic terms for the axions. It is 
possible to retain good accuracy by simply treating the 
Kahler metric as diagonal, in which case the axion decay 
constants can be extracted from the diagonal components 
of the Kahler metric as: 



2Kii , fg 



2K„ 



(B6) 



The decay constants in ( B4 ) arc in fairly good agreement 



with analytical estimate above. 

The next step is to determine the unitary transforma- 



tion lA that diagonalizes the mass matrix ( 14 1. It is given 
by: 



0.706 0.708 -0.019 
UKi [ 0.706 -0.702 0.093 
-0.053 0.079 0.995 



(B7) 



It is convenient to express all the masses relative to the 
gravitino mass scale TO3/2- The axion masses obtained 
from diagonalizing the mass matrix ( 14 1 are: 



rn^i 286 7713 / 2 , 777^2^6.3x10 ^^7713/2, (B8) 
771^3 sa 4.0 X 10^^17773/2. 

Expressing the above masses in units of Planck mass we 
obtain: 



777i « 1.1 X IQ-^^ ml , « 5.2 x 10"^"^ 777^ . 



ml, « 2.1 X 10""^ 7772 , 



(B9) 



7i « 284.9 7773/2 , 7772 ~ 2.0 7r73/2, 7ri3 sa 1.1 7)7,3/2. (B2) where we used the value of the gravitino mass specific to 

the above numerical examplcr] 7^3/2 « 277 GeV. We can 

the 



If we define the initial axion basis (t, &) = {ti,t2,9i), then 
the axion kinetic terms are diagonalized by the unitary 
transformation U: 



now compare the masses of tne axions obtained numeri 
cally with the approximate analytical formula (AS I. 



1.00 -10-4 0.01 
10-4 1.00 0.02 
-0.01 -0.02 1.00 



(B3) 



The above matrix indicates that there is very little mix- 
ing among the components, which agrees with the remark 



^ Such a low gravitino mass scale is an artifact of the toy model 
that has only two moduli. In this case the seven dimensional 
volume Vx = s^^^s"^^^ ~ 3880 is rather large, which makes the 
gravitino mass smaller than in the more realistic examples where 
Vx ~ 0(500 - 1000). 
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Using the general formula (A8), and using the expres- 
sion for /i+2 in (B6) with Ci+2 



1, we obtain an order 



of magnitude estimate for the axion masses 
TO^j 68m3/2 , rh^^ ~ 3 x 10"'^^ 



m^3 « 6 X 10 ^°m3/2 



(BIO) 

which upon comparing with the exact numerical result 
( B8 1 , confirms that ( AS ) is a valid approximation for 



masses of both light and heavy mass eigenstates. 

Let us express the original axion fields (t, 9) in terms 
of the canonically normalized mass eigenstates (before 
taking QCD effects into account): 



ti 

t2 



23.3i/;i 
11.71/^1 



23.4i/;2 

11.6-02 



(Bll) 



6*1 = -O.eV'i + 6.3 X 10"^V2 + 0.8^-3 . 



The effective decay constants fx for the mass eigenstates 
ipK can be computed using the general formula (17), 
where index K now runs over 1 < K < N -\- 1. Keeping 
the coefficients NJ'^^ arbitrary, the effective decay con- 
stants for our toy example are given by: 



Note that according to the arguments in section III the 
mass of the heavy axion eigenstate '0i is the same as that 
for ipi from ( B9 ) while the mass of lightest eigenstate V'a 



in (B9) is completely modified, receiving a mass predom- 



inantly from QCD instanton effects. The mass of -02 also 
receives a noticeable modification from QCD instantons. 
The eigenstates 02 and '03 are therefore modified from 
"02 and 03 respectively. -03 can be identified with the 
QCD axion, which agrees with the expectation from the 
arguments in section |III| 

The eigenstates with significant couplings to the visible 
sector are the heavy state -01 , and the light eigenstate 03 
that picks up its mass from the QCD instanton effects 
and is of m'}^'^ = ©(A^^^/Ja). In more general cases, 
there are generically other eigenstates with masses heav- 
ier than m'^'^^ so they will also couple appreciably to the 
visible sector, in particular to E ■ B. On th e othe r hand, 
eigenstates much lighter than ml'^'^ (-02 in (B14) above) 



will not couple appreciably to E ■ B, as their couplings 
to the visible sector are expected to b e suppressed by the 



mass ratio (to 



qcd\2 



From (B14| we see that this is 



indeed the case for i/)2. 



/i = (l46.5 X Nl' 



h 
h 



146.9 X Nl' 
(3.8 X TV™ H 



+ 73.2 X A^2™ 



9.8 X A^J'" 



(B12) 



TO. 



where we took into account the factor of 27r multiplying 
the fields ti. From the arguments in section |IIH it is 
clear that the lightest eigenstates 02 and 03 with w 



-101 



TO, 



5.2 X 10- 



much smaller than ' 



"i, and ml, « 2.1 x 10 



133 



TO 



p;, which are 



^exp 



10"*2to2; (see eqn.(|2l|)), and 



with effective decay constants fi ^ ©(Mqut), are good 
candidates for the QCD axion. As explained in section 
|III[ the QCD axion is in general a linear combination 
of t,'s, for simphcity we choose Nl'"" = TVJ'" = 1. This 
linear combination representing the QCD axion can be 
written in terms of the mass eigenstates, before the QCD 
instanton effects are taken into account: 



^ qcd 



2tt (A™ti 



A2'"'t2 



(B13) 



27r(ti + t2) ~ 219.8^/'i + 74.10^2 + 5.90^3 : 



Finally, we take into account QCD instanton effects 
and determine precisely which axion candidate satisfies 
all criteria for being the QCD axion. The linear combi- 
nation of final mass eigenstates, which directly couples 
to the visible sector is given by: 

(B14) 



^qcd 



27r(Ar™ii + Nrt2) = 2^{ti + t2) 
w 219.8 0ii -1- 5.5 X 10^28^^ _ 74 3 ^, 

and the final mass spectrum is given by: 

1.1 X lO^^W , TO^- « 3.3 X 10"^°^TO 

^2 ~ 5.5 X 10 



m^7 



TO 



^3 



-''ml 



(B15) 



Appendix C: Axion Decay Constants 

In this section, we compute the generic spectrum of 
axion decay constants fa, as well as effective decay con- 
stants fa that are relevant in determining the axion cou- 
plings to the visible sector, in a framework in which all 
moduli and axions are stabilized from a set of microscopic 
"data" . Instead of relying on any specific choice of the 
G2 manifold, we will use some of the generic properties of 
the Kahler metric for the axions and the moduli to give 
an order of magnitude estimate. We perform a simple 
statistical analysis of the axion decay constants by con- 
sidering a general class of Kahler potentials consistent 
with G2 holonomy where the seven-dimensional volume 
is given by 



M 



N 



V = 47rV3^c,n^^ ' ^af = 7/3, Vfc , (CI) 



k=l 



i=l 



where ~ C'(l) are integer coefficients and the expo- 
nents \ < <\ are multiples of ^ so that each product 
term in the sum contains a maximum of seven distinct 
factors. The latter condition was motivated by the form 
of the Kahler potential for Joyce orbifolds and imposed 
for simplicity but may be relaxed by considering smaller 
values of a\. We assume that a manifold has A^ = 50 
moduli and that the number of distinct terms inside the 
sum is M ~ 0(1000). To determine the moduli vevs at 
the minimum we first need to solve the system of equa- 
tions to determine the parameters explained in [17] . 
Here it is also assumed for simplicity that the integers Ni 
inside the gauge kinetic function of the dominant gaug- 
ino condensates are random sets containing I's and 2's. 
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With these inputs, the distribution of is presented in 
Figure [C] which is rather broad. Such broadening can 
be attributed to the variation of the integers Ni since 
increasing a particular Ni leads to a slight decrease in 
the value of the corresponding observed from nu- 

merical simulations. The distribution of a.i in Figure [C] 
was generated by solving a system of equations \17\ for 
200 distinct randomly generated Kahler potentials of the 
form described above with = 50 moduli. For each 
solution we have verified that all 50 parameters Hi al- 
ways add up to 7/3, as expected. The mean and thermal 
deviation are given by: 



7 
150 



0.047, S.D.(a,) w 0.011. 



(C2) 



600 



400 



200 



0.03 0.04 0.05 0.06 0.07 0.08 0.09 



FIG. 5. Distribution of Ui obtained for 200 randomly gen- 
erated Kahler potentials consistent with G2 holonomy with 
= 50 moduli. The for each case, the integer parameters 
Ni were randomly generated sets containing Is and 2s. The 
mean value is at — ~ 0.047 and the standard deviation 
S.D.(aO 0.011. 



We now use the same 200 sets of parameters Hi along 
with the corresponding integers Ni to compute the mod- 
uli vevs at the minimum by following the analysis in |17j , 
where we choose Q = 60 and Q — P — 3, and then use 
those to evaluate the Kahler metric at the minimum. The 
Kahler metric Kab is diagonalized for each set with all 
eigenvalues strictly positive. Note that we have not in- 
cluded the meson Kahler metric into our analysis. The 
corresponding decay constant is an order of magnutude 
larger that the axion decay costants corresponding to tiS, 
hence its inverse gives the smallest contribution when 
computing the effective decay constants as will be seen 
later. We then compute the axion decay constants by 
using the definition = nip. Two sharp peaks 

colored in red in Figure [C] represent the distribution of 
the axion decay constants on a logarithmic scale. The 
double-peak shape can be traced back to the fact that 
there are two distinct sets of moduli vevs. The peaks 



correspond to the following values of the decay constants: 

(/L)pcak, « 1-3 X 10"GeV, (/L)p,,k. « 2.9 x lO^GeV 
The above result can also be obtained by a simple approx- 
imation. Neglecting the difference due to the integers Ni 
and using the diagonal components of the Kahler met- 
ric Kab = of ' ^ = results in the following 
parametric dependence of the axion decay constants on 
N and Q : 



h 



1.6 



N 



Q 



X IQi^^GeV, 



(C3) 



which for Q = 60 and iV = 50 results in w 1.8 x 
10^^ GeV, very close to the values obtained by the sta- 
tistical method above. It is clear that by keeping the 
number of moduli N fixed while increasing the dual Cox- 
eter number Q of the gauge group, one can lower the 
axion decay constants. Thus, it seems reasonable to ob- 
tain axion decay constants of magnitude 10^^-10^^ GeV, 
consistent with standard gauge unification. 




FIG. 6. The broad distribution in light blue color corresponds 
to absolute values of the effective decay constants on the log- 
arithmic scale Logj^od/Ll/GeV) whereas the double-peaked 
distribution in red color corresponds to the original decay 
constants LogjQ(/i,/GeV) on the same scale. 



Finally, we use eqn. (17 1 to compute the effective de- 



cay constants, where for each set of /lS, the unitary ma- 
trix Ukl used in (17) was obtained by diagonalizing a 



randomly generated symmetric matrix, while Ukl was 
the actual unitary transformation that diagonalized the 
Kahler metric. Here it is assumed that the integers N^^^ 
of the visible sector gauge kinetic function are randomly 
generated sets containing 0, 1, 2. The distribution of ab- 
solute values of the effective decay constants on a loga- 
rithmic scale is presented in Figure [C| in light blue color. 
Since the distribution is clearly non-Gaussian, the peak 
value is somewhat smaller than the mean value: 



(i/li) 



1.6 X 10^^ GeV 



Peak 



(i/li) 

\ / 1 



3.6 X 10^^ GeV 



Thus, from the above analysis the value of the effective 
decay constants are expected to be few x 10^^ GeV, the 
same as the scale of standard gauge unification. 
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